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^ ' Abstract. We perform a high-temperature expansion of scalar quantum field theory on 

fuzzy (DP" to third order in the inverse temperature. Using group theoretical methods, 
we rewrite the result as a multitrace matrix model. The partition function of this matrix 
model is evaluated via the saddle point method and the phase diagram is analyzed for 
various n. Our results confirm the findings of a previous numerical study of this phase 
diagram for <CP^. 
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■ 1 Introduction 

00 ■ 

I Fuzzy spaces are noncommutative geometries which arise from quantizing certain compact 

• ' Kahler manifolds. The most prominent such space is the fuzzy sphere, which was first con- 

. structed by Berezin [1]. In the original construction, the aim was the same as that of geometric 

I quantization, i.e. to provide a general quantization prescription for a particle whose phase space 

is an arbitrary Poisson manifold. Today, fuzzy spaces attract most interest for different reasons: 
First, fuzzy spaces appear quite naturally in various contexts in string theory where they re- 
^ I place parts of the classical geometry of the target space with an approximate quantum geometry. 

^ ' Closely related is the observation that fuzzy geometries seem to emerge from the dynamics of 

^ - matrix models and thus they could be crucial in background independent formulations of theo- 

ries of gravity. And finally one can regulate quantum field theories on Kahler manifolds by 
putting the theory on the corresponding Berezin-quantized or fuzzy manifold. 

The idea of using fuzzy spaces as regulators for quantum field theories goes back to the 
early 1990's OS]. This approach is very appealing, as the definition of scalar quantum field 
theories on fuzzy spaces is under complete control: All functional integrals are automatically 
well-defined because the algebra of functions on a fuzzy space is finite dimensional. Taking the 
large volume limit of the fuzzy space, we can even regulate scalar quantum field theories on 
flat spaces and thus try to compete with the lattice approach. The main advantage of fuzzy 
regularization over the latter is that all the isometries of the original Kahler manifold survive 
the quantization procedure. 

Particularly nice spaces to use in a fuzzy regularization are the complex projective spaces, 
as they are the Berezin-quantizable manifolds with the largest possible symmetry groups. Fur- 
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thermore, their quantization is straightforward and can be done completely in terms of group 
theory. As usual in a "good" quantization, real functions are mapped to hermitian operators on 
a Hilbert space, which is finite dimensional in Berezin quantization. Real scalar field theories 
are therefore simply hermitian matrix models. 

The most prominent hermitian matrix models are given by a potential consisting of a trace 
over a polynomial in the matrix variable. One can therefore switch directly to an eigenvalue 
formulation. In the case of scalar field theories on fuzzy spaces, this is not possible because the 
kinetic term yields a coupling to a number of fixed "external" matrices. 

A first attempt at gaining an analytical handle on fuzzy scalar field theories was made in [4] . 
A new method to overcome the problem of external matrices was then proposed in [5]. Here, 
a high-temperature expansion of the kinetic term in the partition function was performed and the 
resulting expressions could be evaluated analytically via group theoretic methods. It was shown 
that the resulting partition function can be rewritten as the partition function of a multitrace 
matrix model. This partition function can then be computed analytically for both finite and 
infinite matrix sizes using, e.g., orthogonal polynomials or the saddle point approximation. 
For CP^, this computation was performed to second order in the inverse temperature /3 in [5]. 
In this paper, we continue this work and generalize the results to third order in (3 and to 
arbitrary CP". 

One of the motivations for this work is to explain the phase diagram for scalar field theory on 
fuzzy CP^ which has been obtained via numerical methods in [S], see also [7] for a more detailed 
study as well as [8] for a review and further numerical results. The numerical results suggest 
that the phase diagram is invariant under a particular multiscaling. We can therefore restrict 
ourselves to the limit of infinite matrix size, in which we can use the saddle point approximation 
to compute the partition function of our model. 

Further reasons to compute the multitrace matrix model to third order in (5 are the possibility 
to use this result in a similar study of scalar field theory on R x CP^ as well as our intent to 
discuss the link to (deformed) integrable hierarchies in future work. 

In the analysis of the phase diagram, we will focus our attention on the three lowest- 
dimensional fuzzy spaces CP^, and CPp^. In the first case, the goal will be to compare the 
resulting phase diagram with the numerically obtained one. The quantum field theory on the 
second space corresponds in the large volume limit to a scalar quantum field theory of (j)^-iyY>e 
on R^. While admittedly it is not clear what the Lagrangian of the field theory on being 
regularized actually is, this presents an example of both a well-defined and renormalizable four- 
dimensional noncommutatively deformed c^^-theory. The theory on CP"^ could be interpreted 
as a regularization of a non-renormalizable field theory, and one might hope for signs of this in 
the matrix model. 

The paper is structured as follows. In Section [2l we review the construction of fuzzy CP" 
and scalar field theory on this noncommutative space. Section [3] describes the high-temperature 
expansion in detail and the results are presented to order . In Section HI we analyze the thus 
obtained multitrace matrix model for the three lowest-dimensional fuzzy CP" and we conclude 
in Section [5l Conventions, rather technical details and helpful intermediate results are given in 
the appendix. 

2 Scalar field theory on fuzzy CP"^ 

The general mathematical framework containing the quantization of complex projective space 
which is referred to as fuzzy CP" in the physics literature is known as Berezin- Toeplitz quanti- 
zation, see e.g. [9] and references therein for a detailed discussion. In the case of CP", there is 
a shortcut to the general constructions of Berezin-Toeplitz quantization which originates from 
the fact that CP" is the coset space U(n + 1)/U(1) x U(n). We will use this group theoretic 
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approach here, as it has the additional advantage of allowing for simple computations of quanti- 
ties like spectra of quadratic Casimirs and their eigenspaces, which we will need for our further 
discussion. 



2.1 Berezin quantization of CP" 

The Hilbert space M'g which we use in Berezin quantizing CP" is the space of global holomorphic 
sections of the line bundle 0(€) over CP" with ^ > 0. As a vector space, is spanned by 
the homogeneous polynomials of degree i in the homogeneous coordinates z^, . . . , z" on CP". 
Recall that CP" = SU(n + l)/S(U(l) x U(n)), and ^ forms a representation of SU(n + l) which 
is given by the totally symmetrized tensor product of i fundamental representations. In terms 
of Dynkin labels, this representation reads as (£, 0, . . . , 0) and has dimension 



Nni ■■= dim(^) = dim{£, 0, . . . , 0) 



jn + i)] 
nW. ' 



We will find it convenient to map the polynomials to elements of the ^-particle Hilbert space 
in the Fock space of n+1 harmonic oscillators with creation and annihilation operators satisfying 
the algebra [oq,, ajj] = 6ai3 and aa\0) = for a, /3 = 0, . . . , n. We thus identify 

^,^span(at^---atjO)). 
The Berezin symbol map ae : End (M) C°°(CP") is defined as 

ae{f)iz) := {z,£\f\z,£), 
where \z,£) are the Perelomov coherent states, 

The quantization map is given by the inverse of a£ on the set 
of quantizable functions. Explicitly, we have 



(jf(End (J^)) C C°°(CP") 



Zai ' ' ' Zcti Zj^-^ • • • Zfj^ 



\2£ 



1 



it 



aLjO)(0|a/3i •••a/3,. 



Furthermore, aj^{l) = 1 and real functions are mapped to hermitian operators in End (J^e)- 

Note that for CP^, the real part of is given by the spherical harmonics with maximal 
angular momentum i. In general, the endomorphisms End (J^) — split into irreducible 
representations of of SU(n + 1) according to 



l(£)n 



©n 



or equivalently, written in terms of Dynkin labels: 

(A0,...,0)C5(£,0,...,0) = (£,0,...,0)©(0,...,0, 



,o(m,0, ... ,0,m). 



The generators of 5u{n + 1) are represented on End (J^e) by the adjoint action of hermitian 
matrices Lj, and we introduce the quadratic Casimir operator according to 



C2f:= [Li,[Li,f]]. 
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The eigenvalues of C2 are positive and given on the irreducible subspace with Dynkin labels 
(m, 0, . . . , 0, m) bj{3 2m{m + n). The degeneracy of each of these eigenspaces is given by 

2 _ ^2 ^ n{2m + n){{m + n - ly.f 
"'"^ '^r^,m-l (m!)2(n!)2 

Because of C2{cr'[^{f)) = a^^{Af), where f £ T,i and A is the Laplace operator on CP**, it is 
justified to identify C2 with the Laplace operator on fuzzy CP"'. 

The matrices Li represent the generators of su(n + 1) and thus satisfy the algebra 

[Li,Lj] =: \fijkLk, 

where the fij^ are the structure constants of su(n + 1). We choose the Lj such that 

tr(Li) = 0, Lf = clI and tr(LiLj) = ^^^"'^ 6ij. 

(n + Ij^ — 1 

In the adjoint representation i? = (1, 0, . . . , 0, 1), they satisfy the Fierz identity 

LfLf = 5''^6^^ ^ 5"^ , ( 1 ) 

from which we conclude that 

iiR{LiLi) = (n + 1)^ - 1 and iiR{LiLj) = 6ij. 
With the above relation, one readily verifies the following identity for the structure constants: 

hjkhji = 2{n + l)Ski. (2) 
Using the overcompleteness relation for the Perelomov coherent states, 

d^i\z,e){z,£\ = vol(CP")l, 

where d/i = ^ is the Liouville measure obtained from the Kahler form oj yielding the Fubini- 
Study metric, one readily deduces a formula for integration: Given a function f £ the 
integral can be written as a trace over the quantized function a7^{f) G End (J^i): 



I 



, , volfCP") , 1,,,, 



2.2 Quantum scalar field theory on (CPp 

As we are interested in matrix models, it is convenient to switch from the label i of our rep- 
resentations to the label Nn/ and drop the subscript. One should, however, keep in mind that 
only for <DP^, there is an I for every value of N. In the following, we will represent elements 
of End (J^i) by hermitian matrices ^ of dimension N x N. 

In the previous section, we collected all the necessary results for writing down a scalar field 
theory on fuzzy CP". Putting everything together, we arrive at the following action functional! 
on End (JTi): 

5[«>] := tr ($C2«> + r $2 + ^ $4^ /^[^.^ ^^.^ $]] + ^ $2 ^ ^ _ 



'^Note that our conventions for C2 differ from [5] by a factor of 2. 

^We implicitly reabsorbed all volume factors by a rescaling of the field $ and the couplings. 
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As we work with hermitian generators Lj, the quadratic Casimir operator C2 has positive eigen- 
values and for r G R and 5 > 0, the action is therefore bounded from below. This, together 
with the finite dimensionality of End (J^), enables us to introduce the well-defined functional 
integral 

Z:= j m e-^^[*] := j d/i,,($) e"''^!*], (4) 

where d/x £)(<!>) is the Dyson measure on the set of hermitian matrices of dimension N x N. 

Recall that we can diagonalize a hermitian matrix <I> according to <I> = QAQ\ where Q G U(A^) 
and A = diag(Ai, . . . , Xn) is the diagonal matrix of eigenvalues of Under this decomposition, 
the Dyson measure splits into an eigenvalue part and an "angular" integration over U(A^): 

N 
i=l 

where dfiui^) is the Haar measur^ and A (A) is the Vandermonde determinant 
A(A) :=det([A^%•)=^(A^-A,■). 

i>j 

In the case of simple hermitian matrix models consisting of traces (and multitraces) over 
polynomials in <I>, the angular integration is trivial, because tr($") = tr(A'"), and reduces to 
a constant volume factor. The remaining integral over the eigenvalues can then be computed by 
standard methods as e.g. the saddle point approximation or orthogonal polynomials. Here, how- 
ever, the kinetic term contains the fixed external matrices Li which obstruct a straightforward 
translation to the eigenvalue picture. 



2.3 The toy models N = n + 1 on CP"" 

In the case N = i.e. when 1 = 1 and End {^J^) forms the adjoint representation of su(n + l), 

the kinetic term of our model ^ can be evaluated explicitly by using the Fierz identity ([I]). We 
find here that 

tii^C2<^) = p^{Ntr{^') - tr($)tr(cl>)), (5) 

where tr(K) stands for the sum over the eigenvalues of C2 on End {Mi). Note that, as necessary, 
the kinetic term vanishes for <^ ~ 1. We will use this class of toy models for consistency checks 
of our computations below. 



3 The high-temperature expansion 

As it does not seem possible to compute the partition function @ analytically, we perform 
a high-temperature expansion as suggested in [5]. That is, we separate out the kinetic term in 
the functional integral and Taylor-expand its exponential, assuming /3 to be small. As /3 is usually 
inversely proportional to the temperature in statistical mechanics models, this expansion is also 
known as a a high-temperature expansion in the literature. For each of the terms appearing in 
this expansion, the integral over the angular part of the Dyson measure can be performed - in 
principle straightforwardly - using group theoretic methods. The results can be rewritten in 
terms of multitrace terms, and, after putting them back into the exponential of the functional 
integral, one ends up with a multitrace matrix model. 

^That is the unique measure on U(7V) which is invariant under left and right group multiphcation and norma- 
lized according to /d/iH(f2) = 1. 
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3.1 Setup of the expansion 

Let us consider our model ^ on fuzzy CP^ with the dimension of the quantum Hilbert space J^e 
being N. The space of quantized functions End (^) is spanned by the generator^] r^, ^ = 
1, . . . , A^^ of u(A^). We start by rewriting the kinetic term of the action in the fohowing way: 

tr($C2$) = tr ($[Li, [Li,<^]]) = tr{Tf,[Li, [Li,T^]]) tr{^Tf,) tr(^>r^) =: K^,u^^,^u. 

Because of C2IN = 0, we have K^^i^^^^^ = Kmn^m^n, m, n = 1, . . . , iV^ — 1. The expansion of 
the kinetic term in the action now reads as 

2 

and we will restrict our attention in the following to the terms up to order 0{f3^). 

We want to perform the integral over the U(A'^) part of the Dyson measure, i.e. to integrate out 
the angular degrees of freedom in <I>. For this, we decompose the hermitian matrix $ according 
to $ = OAOl", where Q G U(A'') and A = diag(Ai, . . . , Aat). The integrals we have to evaluate at 
order 0{f3^) are thus of the form 

A := f dfiHmYlK^,n,tv{nAQ^T^JtT{nAn^Tn,), (6) 

1=1 

where the essential part in index notation is given by 

I dfinm . . . ■ ■ ■ n^^,^^. (7) 

Various algorithms have been proposed in the literature to compute integrals of the type d?]), cf. 
e.g. [ini [HI [12] • The most involved integral of the form d?]) which we are interested in is the one 
for A; = 3, which is already very difficult to handle by the suggested methods. Fortunately, the 
integrals ([6]) allow for a further simplification [5], which is then accessible via group theoretic 
methods. Using tr(A)tr(S) = tr(A (g) B) and AB (g) CD = {A C){B (g D), we rewrite dSD 
according to 



'^k = j dfini^) Kmrni ■■■K„ 

X tr ( (i7 (g) • • • (g) i7)(A (g) • • • (g) A)(r2"l' (g) • • • (g r2"'')(Tmi 



The idea presented in [5] is now to use the orthogonality relation of the Haar measure (|12p to 
evaluate these integrals. We thus have 

J^k = Kmmi ■ ■ ■ Km^nk ^ ^.^^^^ trp(A (g • • • (g A) trp(rmi (g r„, (g • • 

where the sum is taken over the irreducible representations contained in the tensor product of 
2k fundamental representations of SU(A^). The traces trp are taken in the representation p, and 
we have 

trp(A®---® A) = Xp(A), 

where xpi-^) denotes the character of A in the representation p. Characters of representations 
of SU(A'^) can easily be calculated using e.g. the formulas in [13]. The remaining challenge is 
therefore to evaluate trp(Tmi g) g) • • • (g Tm^ g) Tn^.). 



*Cf. Appendix|X]for our Lie algebra conventions. 
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3.2 The restricted traces trp(-) 

Consider of su(A^), and denote their matrix components by Tm^ , a, (5 

1, . . . ,N. The full trace over the tensor products of matrices in index notation is given by 



tr(Tmi 



J 'm.i 'm2k "oiPl ^a2kP2k- 



To evaluate the restricted traces trp(-), we need to project onto the irreducible representations 

which we do using projectors V)^^ constructed from Young symmetrizers. The technical details 
of the construction of these projectors are given in Appendix [Bj Explicitly, we let the projec- 
tor "Pgj, act onto the indices (3 appearing in the Kronecker deltas to restrict to a representa- 
tion p^*'-'); 

tr,(..) (r^, • • • ® Tm2k ) = r'i^t ' ' ' rZt'^2k^ ' ' ' ^"2./?.. • 

The completeness relation (jl4p for the projectors v'^^'^ translates into the following completeness 
relation for the restricted traces: 



' W) =tr(T-mi ®---®Tm2k), 



which can serve as a first consistency check of the correctness of the calculated projectors V.^^ ■ 
A second test is to verify that each individual restricted trace indeed reduces to the character if 
all the Tm are equal: 

trp(i.,)(A^---® A) = Xp(.,.)(A). 

Let us now compute the combined sums of the restricted traces for each type of Young tableaux 
and contract the Tm with the Kmn to simplify the results. That is, we compute the following 
expressions: 



Km-im2 ' ' ' -^"i2fe— i™2fe trp(i,j) {Tui-i ® '^m2 



m2k-i 



"m.2kJ 



For k = 1 and k = 2, corresponding to the contributions at orders 0{f3) and ©(/J^), these sums 
have already been calculated in [5]. They are 



rrm 



+ 



+ 



20tr(K) - (4 + iV)tr(K)2 


- 2iVtr(/^2) 


24N 

(tr/i:)2 + 2tr 
6 

20tr(K) - (4 + TV) tr(K)2 




- 2Ntr{K^) 


8N 

20tr(K) + (Af-4)tr(K)2 


+ 2iVtr(if2) 


20tr(A') + (Af-4)tr(K)2 


+ 2iVtr(i^2) 


24N 



The lengthy result for fc = 3 is given in Appendix [Cl 
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In calculating these results, we used many identities which we will briefly comment on now: 
First of all, using the Fierz identity for the generators of u(A^) as well as the relations for the Lj, 
we compute that for arbitrary A,B £ u(A^), 

K^^ty:{T^A)tT{T^B) = 2cl tT{AB) - 2ti{LiAL,B), 

K^^triT^Ar^B) = 2cl tr(A) tr(B) - 2tr(Li^) tr(L,5). (8) 
Applying these relations to tr(A') := K^^tr{T^T^) yields 

CL = -^ and dg := (n + lf -1- ^ ^ 



The identities ([8]) allow us to successively rewrite expressions involving K^i, in terms of traces 
over products of the Li, which in turn can be reduced using Lf = clIn and the identity for the 
structure constants ([2]). Some useful intermediate results are collected in Appendix O 

3.3 The multitrace matrix model 

Combining the reduced traces with the characters in the various representations, we arrive at 
the following expressions for the J^fci 

tr(K) ,,2n tr(K) ,,,0 

lOtr(K) (-2(l + iV2) +tr(K)) +4(3-2iV2)tr(i^2) . 

=^2 7 tr(A 

N (-36 + N"^ {-7 + N'^y 



H i p-— ^ — — r i ^ tr(A2) tr(A) 



iV2 |^_36 + iV2(-7 + iV2) 

^ 20{-3+2N^)tr{K) + {30-UN^+N^)tT{Kf+2{l8-6N^+N^)tr{K^) ^ ^^^^^ 

m (-36 + iV2(_7 + ^r2)'^^ 
2(l00tr(if) + (-14 + A^2^^^(^)2^2(6 + iV2)tr 



N [-3Q + N'^ {-7 + N"^)' 



■tr(A2)tr(A)^ 



H — -. — -— tr(A)^ 

Af2 (-36 + 7V2(_7 + ^2)2 



The result for J^s is lengthy and because it can be easily calculated from the list of restricted 
traces given in Appendix [Cl we refrain from presenting it here. Note that all the above integrals 
pass the first consistency check: We have = if A ~ l^r . 

To rephrase the perturbative expansion in terms of an effective action, we re-exponentiate 
the terms. That is, we write 

e-/3(5i+52+53) = 1 _ + _ + 0(^4)^ 

2 6 

where we demand that the Si are polynomials in the eigenvalues of order 2i. As the same holds 
by definition for the we can match both sides order by order and arrive at 

Sl = ^l, S2=^-{Jl-^2). 53 = ^(2^f -3^i^2 + ^3). 

2 D 
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We can now perform the second consistency check of our result and compare the re-exponen- 
tiated action with the toy model N = n + 1 from Section [2.31 In the representation N = n + 1, 
we have 

tr(Er) = 2n(l +n){2 + n), tr(K2) = 4n{l + nf{2 + n), tT{K^) = 8n(l + nf{2 + n). 

Plugging this into the expressions for Si, S2 and 53, we find that Si indeed reduces to the kinetic 
term of the toy model ([5]). Since the terms S2 and S3 vanish as required, our results pass this 
consistency check as well. 

By re-inserting the integration over the Haar measure, we can return from the eigenvalues 
to the full hermitian matrices We thus obtain the multitrace matrix model with action 
S = Si + S2 + S3 with 

trm trm ,_^,2 

etc. The involved expressions for 52 and 6*3 are again lengthy but easily calculated from the 
results given above. 

Altogether, we obtained a multitrace matrix model whose partition function approximates 
the partition function of fuzzy scalar field theory on complex projective space up to order 0{f3^). 
This approximation should be valid in particular for large values of the couplings r and g. 



4 Large AT solutions of the model 

The partition function of the multitrace matrix model we obtained in the previous section can 
now be evaluated analytically for finite using the methods of orthogonal polynomials. As we 
are mainly interested in the phase diagram, we consider instead the large A^ limit and use the 
saddle point method to determine the partition function here. We try to be self-contained and 
present the involved steps in detail. 

4.1 The large N limit 

The phase diagram determined numerically in [6] is invariant under the multiscaling limit where 
A 00 and N'^^g as weh as A^^/^/S r are kept fixed. This justifies to solve our model in the 
large A limit to compare it with the phase diagram. In this limit, the discrete set of eigenvalues 
goes over into a continuous function: We rescale Aj — X{i/N) =: A(x), with < x < 1. The 
traces turn correspondingly into integrals: tr($-') = "^^Xj A dx X{xy . 

The formulas for general n turn out to be very lengthy and difficult to handle. Therefore we 
will restrict our attention in the following to the three projective spaces <DP^ , <DP^ and <DP^. The 
first case n = 1 is interesting, as we would like to compare the resulting phase diagram to the one 
numerically obtained in [6]. The second case n = 2 is a well-defined, four-dimensional quantum 
field theory with quartic potential. The third case n = 3 is interesting as the corresponding 
scalar field theory on is not renormalizable. 

The eigenvalues of the quadratic Casimir, the degeneracy of the corresponding eigenspaces 
and the dimension of the representation (£, 0, . . . , 0) (8) (^, 0, . . . , 0) for the cases n = 1,2,3 are 
listed in the following table: 







(CP^ 


€P^ 


eigenvalues of C2 
degeneracy of eigenspaces 

Nn,i 


2i{i + 1) 
1 + 2^ 
i + l 


2£{e + 2) 
(1 + ^)=^ 
i(£ +!)(£ + 2) 


2i{£ + 3) 
^(l+^)2(2 + £)2(3 + 2£) 
i(^+l)(^ + 2)(£ + 3) 
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As the function Nn/ is only surjective for n = 1, and thus we cannot find an ^ for every value 
of N ^ we will rewrite the multitrace matrix model in terms of £. 

From the table above, we easily evaluate the various traces over K appearing in the action 
of the multitrace matrix model. We have for C-P^: 



tr(J^) = l{l + lY{2 + £), tr(K2) = ^^(1 + if [2 + .^^^ 

3^ 

for CP^: 



tr(i^) = U{1 + if {2 + £)2(3 + I), tr(A'2) = i£2(;^ ^ ^^2^2 + ^)2(3 + 



3 

5 

and for GP^: 



ii{K^) = k^(l + if {2 + if {2, + ^2)(4^2 ^ ^21 - 1), 



tr(ir) = ^£(1 + £)2(2 + ^)2(3 + If {A + £), 
tr(ir2) = + ^)2(2 + ^)2(3 + ^)2(4 + 

tr(J^3) ^ ^^2^^ ^ ^^2^2 + ^)2(3 + ^)2(4 + ^)2(5£2 ^ 20^ - 1). 

In the limit £ — )• oo, the expressions for the various matrix models simplify. Switching to the 
eigenvalue description and using the moments 

dx A"(2;), 

we can write them down explicitly. On the three fuzzy CP"s, we have the models 

^5(n=l) ^ ^^3^^^ _ ^2^ _ ^2^ (^2 _ ^^)2 _ (2c? - ^C,C, + C,f 

+ /3£rc2 + Plgci -fj dx dy log | A(x) - A(y)| , 
/3S("=2)=/3^(c2-c?)-/32^(c?-C2)^ 

- (12c? - 36cfc2 + 2144 + 8ci + 10ci(2c? - 3c2)c3 + 5ci) 

f ^2 £4 |. 

+ l3—rc2 + /3—gc4,-— / dxdy log |A(2;) - A(y)|, 

£3 

- /33_ (2c? - 6cfc2 - 3cl4 + lOci + 6ci(2cf - 3c2)c3 + 3ci) 

Zo 

£3 £3 £6 r 

+ rc2 + /3— 5C4 - — / dxdy log |A(x) - A(y)|, 

where the repulsive log-term arises as usual from exponentiating the Vandermonde determinant. 
The ^-dependence of the log-term is due to the factor A'^^, which in turn originated from 
rewriting the double sum as a double integral. In the above expressions, subleading terms in I 
have been suppressed in each summand. 

As a next step, we have to find the appropriate multi-scaling behavior of the constants /3, r, g 
and the continuous eigenvalues A(x). The coefficients of the log-terms determines the desired 
scaling behavior of the total action. We fix the remaining scalings by demanding that the whole 
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action scales homogeneously and that /3g scales with N'^, as for an ordinary hermitian matrix 
model. We thus find the following rescalings: 







A(x)- 


■>rU(x), 


r —7- i'^r, g - 








X{x)- 


■> A(x), 


r — ;> i'^r, g - 








\{x)- 




r — )• I'^r, g - 





Note that the scalings for CP^ indeed agree with the ones numerically determined in [6] as well 
as the ones calculated in [5]- 

As a final simplification, we note that our theory is invariant under ^ — t- — as the potential 
is even. We expect the eigenvalues to respect this symmetr}{^, and therefore we put all the odd 
moments C2n+i, n G N to zero. Moreover, we replace the integral over x by an integral over the 
eigenvalue density p{X) := We thus eventually arrive at the following three models, which 
we wish to solve: 



dAd^p(A)log|A-/i|p(^), 



(n=3) 



/3 



S/3 
27/3 



256/32 2 



-C2 



405 

72/32 ^ 



+ 2r C2 + 2/3gc4 



C2 + 6r C2 + QPgcA 



dA d/zp(A)log|A-/i|p(/i), 
dA d/U p{X) log I A — /i| p(/i). 



4.2 Solving the models 

We will now calculate the partition functions of our models using the saddle point method: In 
the large (. limit, the path integral localizes on classical solutions, or saddle points, of the actior|2|. 
These solutions, which are valid only for a restricted range of the coupling constants, can be 
easily obtained using standard methods in random matrix theory. We start from the actior|l| 



S[p{X)] = [dX p{X)ViX) - [ dA dp p{X) log I A - p\ p{p) + i( I dX 
Jx Jxxx \Jx 



p{X) - 1 



where I is the union of open intervals on the real line over which p{X) has support. The saddle 
point equation is obtained by varying the above equation with respect to p(A): 

ViX)-2j^dpp{p)log\X-p\+C = 0. (9) 

Note that our potentials satisfy V{X) = at A = and we can therefore determine the Lagrange 
multiplier ^ by solving the saddle point equation at this special point if G X: 



Coex = 2 j^dpp{p)\og\p\ 



otherwise, one has to choose a different value of A to obtain ^. We define the free energy F as 
F := — \og{Z), where Z is the partition function of our model. In the saddle point approxima- 
tion, this reduces to F = l3S[p{X)], which we can evaluate using 

F = \ j^dXp{X)V{X)-\i. 



^As we will see later, this assumption is not correct in a small part of the configuration space. At this point, 
it serves as a very useful approximation to keep the terms in the action manageable. 

®Note that we had to switch to the eigenvalue formulation of the matrix models first, as the zero modes 
corresponding to the angular degrees of freedom contained in "3? would have rendered the approximation invalid. 

^We have included a Lagrange multiplier ^ to fix the normalization of the eigenvalue density. 
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To find the eigenvalue density p{X), it is convenient to replace ([9]) with its derivativ^ with 
respect to A: 

nA) = 2/d^^. (10) 

This is a singular integral equation, and its general solution can be found e.g. in [T3], see also |15j . 
First of all, one introduces the resolvent W{X), which is an analytic function on C\X, defined 
according to 

W{X) := 



X — fi 



Note that for large A, we have W{X) ~ j. The resolvent is related to the eigenvalue densi- 
ty p{X) and the Cauchy principal value appearing in the equation of motion through the Plemelj 
formula, and we arrive at 

V'{X) = W{X + ie) + W{X - ie). 

The first equation determines p{X) in terms of the resolvent, and the second equation is a much 
simpler equation than (jlOp . which fixes the resolven10 and thus the eigenvalue density. One can 
show that the resolvent satisfies the Schwinger-Dyson equation 

W\X) - V'{X)W{X) + iRiX) = 0, 

where 

y'{x)-v'{p) 



R{X) = A I dfipip)- 



p — X 

is a polynomial of degree d — 2. The solution to the above equation reads as 



' V ' 

:=a;(A) 

where uj{X) describes the part of W{X) containing the branch cuts. 

Explicit solutions are now obtained by making assumptions about the support I of the 
eigenvalue density p{X). The simplest assumption is that I consists of a single interval. This 
is expected if the potential either consists of one deep well or if the eigenvalue filling is such 
that all the local minima of the potential are more than filled up. In this case, the resolvent 
has to have a branch cut over X := {6i, 62) and the corresponding solution is therefore known as 
a single-cut solution. The resolvent's singular part has to contain exactly two roots, cf. e.g. |16j : 

uj^{X) = m2(A)(A - <5i)(A - 62) = V'^{X) - R{X). 

One can now make a general ansatz for the polynomials M(A) and R{X). Together with the 
self-consistency condition that all the moments c„ satisfy their defining relation 

dAp(A)A", 



When doing this, one obviously has to vary each moment: Sc2 = 2c25c2, etc. 

^Strictly speaking, it fixes the resolvent only up to regular terms, which, however, are absent as can be seen 
from the large A behavior W{X) ~ j. 
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we can solve for all unknowns and determine p{X). Note that the normalization condition on 
the eigenvalue density cq = 1 is equivalent to the less involved condition that the asymptotic 
behavior of the resolvent is W{X) = j + 0{jj). 

When having a double well potential, we also expect solutions where X is given by the union 
of two disjoint intervals I = ((5i, (52)U (ei, £2)- Correspondingly, the singular part of the potential 
contains four roots and we make the ansatz: 

uj\X) = M\X){\ - 5i){X - 52){X - ei)(A - £2) = V'\X) - R{X). 

This solution is known double-cut solution. 

It is important to stress that in general, all solutions will be valid only on a subset of the full 
parameter space of the model under consideration. This subset is characterized by the condition 
p{X) > (and therefore C2n > 0) as well as the condition that X is of the assumed form. It is 
called the existence domain of a solution, and its boundary in parameter space can correspond 
to a phase transition. 

In the following, we will present all the solutions for the various models together with their 
existence domains. For <DP^, we also give the explicit expressions for the free energies. We will 
consider three kinds of solutions: the symmetric single-cut solution, the symmetric double-cut 
solution and the asymmetric single-cut solution. In the latter case, we should strictly speaking 
include all the odd moments C2n~i, n g¥\, which we dropped in our actions. However, the full 
action would be very difficult to handle analytically, and we hope to make at least qualitative 
statements with our truncation. 

The solutions for closely related models, in which r is kept fixed, the coefficient of is 
a parameter and the coefficient of vanishes, have been computed in [17j for the symmetric 
single-cut type and [H] for the two other types. 

4.3 Solutions of the model on CP^ 

For the symmetric single-cut case, we assume that the eigenvalue density p{X) has support on 
the interval X = (— d, +d). The solution we obtain from the procedure described above together 
with the conditions W{z) ~ ^ and Jj-dXX^ p{X) = C2 reads as 



3{2d^rP + 3d^gP + 2d^/3-4:) ^/ - X^{4 - d^ gP{d^ - 4X^) 

= w ' ^ ^ ^ ' 

48 + d^/3{d^g/3^ + 4d^{l3 - 9g) - 24(1 + r)) = 0. (11) 
The free energy of this solution is given by 

= ^ (40 + d^(3 {-12d^g + 4(1 + r) + A(l + r)f3) - 64 log (^0 ) . 

The solution exists, if both p{X) and C2 are nowhere negative. The condition p{X) > amounts 
to 

-, 2(/3 - 35) 



3^ 

while C2 > is always satisfied in the existence domain of the solution (|lip . 

Next, we assume a symmetric double-cut support for p{X) on X = (— t/s + d, —y/s — d) U 
(Vs — d, s + d). The solution here reads as 

p(A) = -g/3AV(^ + d-A2)(A2-s + d), C2 = d=-^, ' = ^TW^V 

vr 2(/3 - 35) 
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To evaluate the free energy, we compute ^ at A = y/s and use the relation 

F = i (j^dAp(A) (y(A) - log(A -V~s)- log(A + ^s))^ + iy(V^) + lV{-V^s). 

The result is 

9g - 3 (3 + 4r + 2r2) /3 + (6g - 2/3) log(4g/3) 
24<7 - 8(3 

The symmetric double-cut solution exists, if s = C2 > 0, i.e. if r < — 1 and 3g > (3 or r > —1 
and 3g < (3 and if s > d. The latter condition yields 

2(/3 - 3^) 

r<-l+ 



and we see that the boundary of the existence domain of the symmetric double-cut solution 
matches that of the symmetric single-cut solution. We therefore expect a phase transition 
at this boundary. At the point (ro,go) = something interesting happens: Here, the 

equation for s becomes trivially satisfied and s is unconstrained. The two cuts can thus be 
arbitrarily far apart. At this point, the action reduces to 

and we have a competition of the single-trace and the multitrace potential term. 

The asymmetric single-cut solution with support on X = (s — d, s + d) with s 7^ is given by 



2gf3y/d^-{s-X)^X{s + X)-d^) 3{l + 3d^g + r + 2gs^) 
P^^^ = n ' = 2-p 



s = y ^gjy ^ 4/3 + g (-12 - 12^2(1 + r)(3 + 5d'g(3' + d^(3{-A5g + 11/3)) = 0. 

To evaluate the free energy, we determine the Lagrange multiplier ,^ at A = s. Our definitions 
then yield 

_ 8(1 + r) + d'^g (-24 + d^/S {-3d^g + 14(1 + r) + 5d'^g{l + r)/3)) + 32^2^, log (|) 
~ 32d'^g ■ 

The asymmetric single-cut solution exits if /? > and if d is real and positive. The first condition 
implies 

-135^/TEg - 135^^^+41^/15/3 
135^V? ' 

and C2 and s are automatically positive. The second condition amounts to g > ^. 
4.4 Phase structure on CP^ 

If existence domains do not overlap, we expect a phase transition at the boundary. If, however, 
two solutions exist for the same parameters, then the solution with the lowest free energy will 
be adopted. The resulting phase diagram is depicted in Fig. [TJ 

Here, the symmetric single-cut, the double-cut and the asymmetric single-cut solutions are 
labelled as I, II and HI. The boundary of the existence domain between I and II describes the 
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r 

Figure 1. The phase diagram on <DP^ for /3 = i. Daslied lines describe phase boundaries, solid lines 
boundaries of existence domains. See the text for more details. 

usual second order phase transition of the hermitian matrix model. The existence domain of III 
is fully contained in the existence domain of II. There is indeed a region of the parameter space, 
where the asymmetric filling yields a lower value for the free energy than the symmetric filling. 
This is particularly interesting, as it was very difficult to extract the II-III phase transition from 
the numerical data in [6]. We can thus confirm the numerical findings. Furthermore, there is 
the forbidden region in the parameter space for g < ^. We expect that higher-order corrections 
in /3 would deform this boundary. 

Altogether, we have obtained the general features of the phase diagram found in [6]: We have 
three distinct phases, which come together at the point (7*0,50) = 1/6), which corresponds 
to (6, c) = (—0.5, 1/12) in the conventions of [6]. This compares to numerically found values of 
(6, c) = (—0.8 lb 0.08, 0.15 ± 0.05). The discrepancy is due to the fact that the triple point is in 
a region of the parameter space, where the kinetic term is not small compared to the potential 
terms. 

The effects due to the asymmetric single-cut region have to be considered only qualitatively, 
because we have dropped all the odd moments from the action using symmetry arguments. This 
was done to keep the solutions under analytical control. The critical line found in [6] corresponds 
here to the dashed curve. Including the odd momenta would presumably straightened this curve. 

The discrepancies compared to [5] arise from the fact that there, a contribution to the action, 
labelled K~'K, was neglected in the large N limit while we included it here. This yielded 
a different model with the opposite sign of the C2 term. 

4.5 Solutions of the model on CP^ 

Let us now be brief in repeating the analysis for CP^: The single-cut solution with support on 
I = (— d, d) is given by 



_ 2/3Vd2 - A2 (-120C2/3 - 1284^^ + 45 (4 + 3£g + 3r + 6gX^)) 
P^^' - 135^^ ' 

C2 = 1(2^2 + d^^/?), -270 + d^(3{405d^g + 270r - 8(-45 + 2c2/3(15 + 16c2/3))) = 0, 

and the boundary of its existence domain is 

4 [2^ 64/3 + 6OV2P 

T 13% ■ 
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The double-cut solution with support T = {—y/s + d, — V s — d) U {y/ s — d, y/ s + d) reads as: 



P(A) 
d = 



2AJ<7/3(2-45/3(s-A2)^) 



vr 



C2 



-1285^ + 45(4 + 3r + 6gs)/3'^ - 120s/3^ = 0, 



and this solution is admissible, if s is real and s > d. This yields the following bounds: 



-4055^ + 3605/3 - 592/3^ 
384/32 



25 64/3 + 60^2^ 



1355 



Note that the left and the right bound touch at a single point in the r-5-plane. 

The asymmetric single-cut solution has support I = {s — d, s + d) and is given by 



Ag^y/d^ - (s - A)2 (A(s + A) - d^) 



1 



C2 = ^d^g{16s'' + lOd^s^ - d 



vr 



'2 + 3^4^/3 



This solution is valid for 



45 (4 + 9d'^g + 3r + 655^) - 120c2/3 - 128c2/3^ = 0. 



4 71% 82^ 26896/3 
^ ~3 ~ 2771% 1822% ■ 

The moment C2 and the center of the cut s are automatically positive. Contrary to the case 
of <DP^, there is no upper bound on r for this solution. 

The expressions for the free energies on <CP^ are not presented as they are lengthy but can 
be calculated quite straightforwardly. 




4 r 



Figure 2. The boundaries of the various existence domains on (DP^ for /3 = i. See the text for more 
details. 

The boundaries of the various solutions are depicted in Fig. [2j The solid line corresponds 
to the usual matrix model phase transition, i.e. the boundary of the existence domain of the 
symmetric single-cut solution. It is also the upper-right boundary for the existence domain of 
the double-cut solution, whose lower-left boundary is the dashed line. The dotted line is the 
lower-left boundary of the existence domain of the asymmetric single-cut solution. The area in 
the lower left corner is the forbidden region in parameter space. 
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4.6 Solutions of the model on (DP 

We list now the three solutions for the model on <DP"^. This model differs from that on CP^ 
only in the magnitude of the coefficients in the action. The supports are chosen in the same 
way as on CP^ and CP^. 

Symmetric single-cut-solution: 



_ 3/3Vd2 - A2 (lOr - 3(-5 + 6c2/3(l + 4c2/3)) + Wg {d^ + 2X'^)) 

PlA) — , 

ovr 

C2 = 1(2^2 + d^g(3), 10 - 3d^(3{15d^g + lOr - 3(-5 + 6c2/3(l + 4c2/3))) = 0, 



-5g 9 + 2J^ +16/3 + 6V6P 
r > . 

Double cut solution: 

2xJgf]{6-mgf3 {s-X'^f) 

pW = , 

vr 

C2 = s, d=^=, lOr + 205s - 3(-5 + 6s/3(l + 4s/3)) = 0, 
2 , 1 - 11 fil /?2 -5<7 f 9 + 2,/f U 16/3 + 6V6P 



-1005^ + I8O5/3 - 1161/3 

< r < 



720/32 ^ ■ ^ 30g 

Asymmetric single-cut solution: 



^ V2 + 9^4^/3 



2d./W9 ' 



30(i2c, + lOr 2{)gs^ - 3(-5 6c2^(l + 4c2/3)) = 0, 



3 Jbg AlJp 1681/3 



2 ^/2;5 10^105 4509 




4 r 



Figure 3. The boundaries of the various existence domains on <CP^ for 13 ~ \- 

The boundaries for the existence domains are presented in Fig. [3l The meaning of the lines 
is the same as for CP2. Not surprisingly, the phase diagram is essentially identical to that 
for CP2. Unfortunately, there is no feature hinting at the non-renormalizability of (/)^-theory 
on 1R6. 
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5 Conclusions 

In this paper, we computed the partition function of scalar quantum field theory on fuzzy CP" 
to third order in the inverse temperature /3, generalizing the results of [5]- As this theory can 
be interpreted as a noncommutative deformation of scalar quantum field theory on R^'^ in the 
large limit, we also demonstrated the existence of a nontrivial such theory. 

We started by expanding the exponential of the kinetic term in the partition function. We 
then used group theoretic methods to integrate out the zero modes of the action and obtained 
a multitrace matrix model. This model was then solved via the saddle point approximation in 
the large limit. In principle, however, the partition function of the matrix model could have 
been computed as well at finite N using orthogonal polynomials. 

We presented the explicit classical solutions on which the partition function localizes in the 
large limit and discussed the arising phase diagrams for CP^ and CP^. We confirmed 
the findings of the numerical analysis of [SJIT] for CP^ and reproduced qualitatively - and partly 
quantitatively - the phase diagram found numerically. That is, we confirmed the existence of 
three distinct phases and confirmed analytically their properties suggested by the numerical 
studies. We also found a triple point which agrees to an acceptable degree with the one found 
numerically. 

Here, it was particularly interesting that we found a large region of the parameter space in 
which an asymmetric single-cut solution was energetically favorable to a symmetric double-cut 
solution, even though our potential was symmetric. Such situations have been studied in the 
past, see e.g. [19j . Physically, the existence of this spontaneous symmetry breaking in the large A'^ 
limit can be explained as follows. Consider the matrix model at finite N. Here, we have to 
introduce explicitly a symmetry breaking term for an asymmetric phase to exist, as tunnelling 
of the eigenvalues would otherwise restore symmetry in the eigenvalue filling. After taking A^ 
to infinity, there is no more tunnelling, and the symmetry breaking term can be safely switched 
off, preserving the asymmetric configuration. 

It would be interesting to push the analysis of the phase diagrams further and, for example, 
to include all odd moments and examine the possibility of a smooth transition of the filling 
fractions between phases II and HI. We then might be able to reproduce the slope for the linear 
phase boundary found in [6]. It might also be interesting to compare our results to the findings 
of [20], where the phase structure of noncommutative field theories on Moyal space was analyzed 
and |21j . where questions arising from [20] were discussed on the fuzzy sphere. Moreover, we 
intend to use the results found here to study scalar field theory on R x 5^ as well as the relation 
of our multitrace matrix model to (deformed) integrable hierarchies in the future. Finally, recall 
that multitrace matrix models had been proposed as candidates for conformal field theories with 
c > 1 coupled to gravity [T7]. One might be able to make sense of our models in this context, 
as well. 



A Lie algebra conventions 



Note that our conventions differ slightly from those of [5]- Everywhere in our discussion, we use 
orthonormal hermitian generators = r^, ;U = 1, . . . , A^^ of u{N), which satisfj0 

tr (r^r.) = ^ and r^^r^^ = 5^' 6^^. 

The generators of u(A^) split into the generators r^, m = 1, . . . , A''^ — 1 of su(A^) and t^2 = 
-y^lj^. For the former, the Fierz identity reads as 



'm ' m u V N 
^We always sum over indices which appear twice in a product, irrespective of their positions. 
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The Haar measure of SU(A^) satisfies the following orthogonality relation: 



kl 



dim(p) 



(12) 



where Q G SU(A^), p is a finite-dimensional, unitary, irreducible representation and denotes 
its complex conjugate, see e.g. [5] for the proof. 

B Projection onto irreducible representations 

Given a Young tableau A with rows Aj describing a partition of d, we define the symmetrizer s\ 
and the antisymmetrizer ax as operations symmetrizing over rows and antisymmetrizing over 
columns: 



sx : V'^'^ 
ax: V^"^ 



• (g) Sym^'^y c 



where V is the fundamental representation and fii are the rows of the conjugate partition to A. 
The Young symmetrizer cx is now defined as 

CA := sxax, 

and forms a projectors c| = acx with a being the product of all hook lengths of the Young 
tableau A. For more details, see e.g. [22] . To obtain all irreducible representations, we have to 
consider all standard Young tableaux, i.e. Young tableaux, in which the numbers in the boxes 
increase both downwards and to the right. At order 0{/3), we encounter the tableaux 



2 




1 


2 



while at order we need to consider 



2 


<8) 


3 


<8) 


4 




1 


2 


3 


4 



e 


1 


2 


CO 




1 


2 


4 




1 


3 


4 




4 




3 




2 





1 


2 


e 


1 


CO 


3 


4 




2 


4 



A complication sets in for decomposing tensor products of more than four fundamental rep- 
resentations into irreducible ones in this manner: Young symmetrizers of Young tableaux which 
correspond to the same Young diagram are no longer mutually orthogonal and have to be or- 
thogonalized by the Gram-Schmidt method. The 76 Young tableau^o which originate from the 
tensor product |T] ■ ■ ■ 



1 : 2 : 

6 



6 correspond to 11 different Young diagrams: 
3 : 



1, 7. ™ 



T2, A-.^TTD, 5: 

rm 
9 : - , 10 



11 : 



Without spelling them out explicitly, we will order the Young tableaux of a given type lexico- 
graphically. That is, a Young tableau precedes another one, if it has a lower number in the first 
box in which they differ (going through the tableaux from left to right and top to bottom). We 



'^Most of the calculations in the following are done using a computer algebra program. 
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then label them by their type and their position in the ordering: The tableau A^^'^\ for example, 
is given by the third tableau of the second type: 



A{2,3) 



1 


2 


5 


3 


4 


6 



(i i) 

The projectors Vq onto the irreducible representations are again given by the Young sym- 
metrizers c^i^j), except for the ones which are not orthogonal. For example, c^{2,i) c^{2,5) / and 
therefore we need to define 



-0(2,1) 1 1 

'6 — C;^(2.i) 2'^A(2.i)C;^(2,5), 

02 Or, 



(13) 



where the are the products of the hook lengths of the Young diagram of type i, such that 
c^^(i,j)C)^(iJ) = aiC^(i,j). The same is the case for the following pairs of Young symmetrizers: 

Ca(3.1)Ca(3,8), C;^(3,2)C;^{3,9), C;^(3,3) C;^{3,9) , C;^{5,1) C_)^(5,5) , C;^(6,l) C;^(6,ll) , 

Ca{6,1)C;^(6,12), C;^(6,2)C;^{6,13), C;^{6,2) C;^{6,14) , C;)^{6,3) C;^(6, 15) , C_;^(6,4) C;^{6,16) , 

Ca(6,5)C_)^(6,15), C;^(6,7)C;^{6,16), C;^{8,1) C;^(8,7) , C;^{8,1) C;)^(8,8) , C;^(8,2) C;^{8,9) . 

After defining the projectors V^^'^^ appropriately as in (jl3p . they are orthonormal: 

Moreover, the projectors now give the decomposition into irreducible representations: If we have 
a tuple of 6 objects (ai, . . . , og) onto which the Young symmetrizers act, we have 



(ai, ... ,06) = yZ'P6'^\^i^ ■ ■ ■ ,"6) 



(14) 



C Detailed results at order 0{f3^) 

The contracted restricted trace for /c = 3, 



read as: 

"n 



rm 



720iV2 tr(i^) V + ^^^^^ + + + ^^^^ + 

- 8N'^ tr {K^f + 2N tv{Kf (-80(1 + nf + (32 + 3N) tr {K^)) 
+ 16N^ tr{K) (4(1 + n) tr (i^^) + (^-3^^ ^ ^ 

" UANtriK) ( " ^^^^ + + + + 

+ 2tr(K)2 (l6(l + n)^ + (8 + 9iV)tr (K^)) + 24iVtr(i^) tr {K^) 



2\2 



SOAf 



4(1 + n) tr(ir)^ + (4 + A) tr(K)'^ 



+ 2tr(A:) ((8 + 3A)tr (A:^) - 16(1 +n)2) +8A^tr (K^ 

^ ^4(40 + 3A)(1 + n) tr(A:)3 - (40 + A^(12 + A)) tr(A)^ 



144 A2 tr(A') 
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+ 8iV^ tr {K'^f + 2N ii{Kf (80(1 + nf - (32 + 3iV) tr [K'^)) 
- leiV^ tr(K) (4(1 + n) tr [K'^) + tr (K^^ 
1 



2\2 



144iV tr(i^) V^^^^ + + ^("^ + ^) " 

- 2 tr(K)2 (16(1 + n)2 + (8 - 9iV) tr {K^)) + 24iV tr(JO tr {K^) ) , 

16 (2tr(K)2 (4(1 + n)2 - tr (K^)) + tr (k2\2^ 



TD 



2\2 



45iVtr(A') 

- 3g^2fa(j^-) (80(1 + ^) + (-20 + tr(i^)' + 47V' tr {K') 

+ 2Af tr(/^)2 (16(1 + nf + (-4 + 3iV) tr {K^)) 

3 



80iV 



4(1 + n) tr{Ky + (-4 + N) tr{Ky 



+ 2tr{K) (16(1 + n)2 + (-8 + 3iV)tr {K^)) +8iVtr {K^ 



T2 



2\2 



80(1 + n) tr{K f + (-20 + N^) tT{Ky - AN^ tr {K^) 



36A^2 tr(/s:) 

+ 2Af tr(/C)2 (-16(1 + nf + (4 + 3iV) tr {K^)) 
+ 4Af2 tr(i^) (-8(1 + n) tr (K^) + tr {K^^ 



4(-40 + 3iV)(l + n) tr{Kf + (40 + (-12 + N)N) tiiKy 



■ 144iV2tr(i^) 

+ 8iV' tr {K'^f + 2NtT{Kf (80(1 + nf + (-32 + 3N) tr {K^)) 
+ 16iV2 (4(1 + n) tr (i^^) + tr {K'^)) ) , 



1 



720A^2 



4(-40 + 3iV)(l + n) tr{Kf + (40 + (-12 + N)N) tr{Kf 



+ 8Af3 tr (/C^^^ _^ 2A^ tr(i^)2 (80(1 + n)^ + (-32 + 3iV) tr {K'^)) 
+ 16iV2 tr(K) (4(1 + n) tr {K'^) + tr (K^^ 



In the following, we present some more details needed in the calculation of the contracted 
traces. We start by computing the following traces over the generators of su(n + 1) in the 
A^-dimensional representation: 



tv^LiLjLiLj) 



ti{Kf 



(n + 1) 



tr(i<:) 



tic{LiLjLi^LiLjLj^) 
ti{LiLjLf^LjLiLk) 



tr(K)(8(l + nfN^ - 6(1 + n)N'^ tv{K) + tv{Kf 



8N^ 

tr(K)(-2(l + n)N^ + tr{K)f 



8N^ 



1 



tiiLiLjLk)ti{LjL,Lk) = --tr{K') 



tr(K)3 3tv{K)tv{K'^ 



4iV4 



+ 



tr{LiLjLk)tr{LiLjLk) = tr(LiLjLfc) tr(LjLjLfc) + 



8iV2 

{l + n){tr{Kf - N^tr{K^)) 
4iV2 • 



22 



C. Samann 



Using these, we obtain 

tr(Er)2 

tiiK) 

K^,juiK^2i^2 tr(r^iT/,2ri,,T,,2) = 2tv{[Li, Lj][Li, Lj]) = -2(n + 1) ^ , 



Ul V2 ^tJ.i 1^3 (tTji Tjyj ''T'S ''"l'2 ''^3 ) 

-^/^l Ui V2 ^tJ.3 1^3 (17*1 Tjyj T^2 ''T'S ''"I'S ''V2 ) 



2(n+ l)tr(K)2 

\2- 



-^/^ii/i-f^/^2!^2-^M3i'3 ^'^(''7*1^/^2 ''"i^i'^Vs ''"1^2 ''V3) — 2(n + 1) ^tr(iir ) ^ 5 

Kt^JUlKf,2^2K^^l3<^3 ^'^{Vt^2Vyi^^2Tu3) = tr(J^3) + 2(n + 1) (^tr(iC2) - ^ 
as well as 

8(1 + n)2tr(K) 



-^/^l V2 ^ii3 ^3 (''7*1 ^A'2 ''7*3 ) (''"i^l ''"i^2 '^Vz ) 

-^/^l !/l -^/^2 ^2 -^A*3 ^^3 ^-"^ (''7*1 ''7*2 ''7*3 ) ^"^^'^Ul ^^^3 '^^2 ) 

^il\Vx^il2f2 -^A*3 ^^3 (''7*1 ''7*2 ) tr (Ti^-^ rjy2 T^3 Tjyg ) 

-f^/xi i/i -f^/x2 V2 -^A*3 1^3 (''7*1 ''7*2 ) tr (t^i T^2 ''"1^3 ) 



2tr(K)tr iVtr(K2)' 



N tr(i^) ' 

tr(K)tr (^2) 

iV ' 
4(1 + n)2tr(K) 



These expressions can be easily checked for n = 1 and = 2, 3, 4 by comparing them to the 
results of a direct computation. 
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